Nonlinear magnetoinductive waves and domain walls in composite metamaterials 
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We describe novel physics of nonlinear magnetoinductive waves in left-handed composite meta- 
materials. We derive the coupled equations for describing the propagation of magnetoinductive 
waves, and show that in the nonlinear regime the magnetic response of a metamaterial may become 
bistable. We analyze modulational instability of different nonlinear states, and also demonstrate 
that nonlinear metamaterials may support the propagation of domain walls (kinks) connecting the 
regions with the positive and negative magnetization. 
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The past decade observed many advances in the de- 
sign and engineering of artificial structures with unique 
electromagnetic response, which have broadened signif- 
icantly the range of possible wave phenomena that can 
be accessed in experiment. In particular, it has been 
shown that the composite structures may allow realiz- 
ing the materials with simultaneously negative dielectric 
permittivity and magnetic permeability, also known as 
left-handed media |l|, the unique materials because of 
their surprising and often counterintuitive electromag- 
netic properties. The composite metallic structures con- 
sisting of arrays of wires and split-ring-resonators (SRRs) 
have been demonstrated to possess left-handed properties 
in the microwave frequency range Q ■ 

A standard theoretical approach for analyzing the 
properties of composite metamaterials is based on the 
effective medium approximation for both linear 0- Q 
and nonlinear 0, Q metamaterials. In such approach a 
microstructured composite is treated as a homogeneous 
isotropic medium characterized by effective macroscopic 
parameters. This approximation is justified when the 
characteristic scale of the wavelength of the electromag- 
netic field is much larger than the period of the mi- 
crostructured medium. Moreover, this effective medium 
approach is based on a simple averaging over the lat- 
tice of micro-elements, and usually it does not take into 
account any internal (or eigen) modes of the structure, 
which might exist due to interaction between the indi- 
vidual elements via microscopic electromagnetic fields. 
Such eigen modes in the material containing arrays of 
SRRs are referred to as magnetoinductive waves @. 

In this paper we study a novel physics of nonlinear 
composite metamaterials stipulated by magnetic interac- 
tions between SRRs. We consider a cubic lattice of par- 
allel SRRs and demonstrate that the effective coupling 
between the resonators is highly anisotropic. We derive 
discrete coupled equations and describe the properties of 
linear and nonlinear magnetoinductive waves. We show 
that in the nonlinear regime the magnetic response of 
a metamaterial may become bistable, and we analyze 
modulational instability of different nonlinear states. We 
demonstrate that nonlinear metamaterials may support 
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FIG. 1: (color online) (a) Three-dimensional geometry and 
(b) side view of a cubic lattice of split-ring resonators. Lines 
in (b) show the direction of microscopic magnetic field. 



the propagation of domain walls (kinks) connecting the 
regions of positive and negative magnetization, and sim- 
ulate their dynamics numerically. 

Model. We consider a three-dimensional cubic lattice 
of identical parallel SRRs, as shown in Fig. [Ja). We 
assume that the SRR slits are infilled with a nonlinear 
dielectric with a Kerr-like nonlinear response, i.e., with 
dielectric permittivity e = e; + a|i?| 2 /|i? c | 2 , where E c is 
the characteristic nonlinear field, and a — ±1 correspond 
to the focusing and defocusing nonlinearity, respectively; 
the case a = describes a linear response. 

An external electromagnetic field induces currents in 
SRRs; this generates additional magnetic fields which 
determine an overall magnetic response of the compos- 
ite structure. In this paper, we take into account the 
interaction between the nearest neighboring magnetic 
resonators induced by microscopic fields, as shown in 
Fig.flpj), This approximation describes the basic physics 
due to the SRR interaction, since the magnetic field of 
each SRR decays fast enough and the interaction between 
the next neighboring magnetic resonators is weak. In a 
linear regime, it is possible to take into account inter- 
actions between all resonators in the metamaterial sam- 
ple However, in the nonlinear problem the nearest 
neighbour approximation simplifies the solution. 

We use the indices n, q, m to denote the position of the 
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resonators along the axis x, y, z, respectively. Each SRR 
represents an effective oscillatory circuit with inductance 
L = Ana [ln(8a/r) — 7/4] / c 2 of the loop, nonlinear capac- 
itance Cnl ~ e (\E g \ 2 )r 2 /Adg of the slit, and resistance 
R = 2na/aS c fi of the SRR wire (see, e.g., Refs. 0,0 and 
references therein) , where a is the radius of SRR, r is the 
radius of the wire, d g is the size of the SRR slit, E g is the 
electric field induced in the SRR slit, a is the conductiv- 
ity of the SRR wire, S c s is the effective cross-section of 
the wire, and c is the speed of light. 

In the system of magnetically interacting SRRs, we 
also have an additional mutual inductance, described 
by the matrix M, acting between the neighboring res- 
onators. The external electromotive force £ n .q,m in each 
oscillatory circuit is determined by the external magnetic 
field H , £ n , q . m = -{ira 2 / c)(dH / dt) n , q , m . Electromo- 
tive force induced in the resonator at the site (n, q,m), 
due to the mutual inductance with the other resonators, 
can be written as follows, 
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where M„ 



is the elements of the mutual in- 



ductance matrix describing the interaction between the 
resonators at the sites (n, q, m) and (n', q' , m'); Mo,o.o 
. Using symmetry, we can present all non-zero elements 
of the matrix M as follows 



-Mi, o,o = Af-1,0,0 = M i.o = M _i.o = Mi 



(2) 



Mo )0 ,i = Mo,o-i = --M||, (3) 

and derive a set of coupled equations governing the dy- 
namics of the electric currents in SRRs 



dt 



n,q,m ~t~ ^n.g.m — &n.q,m 
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where U n ,q, m is the voltage across a slit of the corre- 
sponding resonator. Capacitance of the SRR slit is a 
sum of linear and nonlinear parts, CNL(|^n.<j.m| 2 ) 
( " + ACNL(|C/n,g,Tn| 2 )j and we assume that the nonlinear 
part is much smaller than the linear one. 

We assume that the fields vary harmonically in time 
[i.e., ~ exp(iwt)] neglecting the generation of higher har- 
monics, and use the approximation of the slowly-varying 
amplitudes for analyzing Eqs. (0} and JSJ). The magnetic 
momentum of each SRR is proportional to the current 



in this resonator, m„ !( j, m = Zo7ra 2 X„. g , m /2c, where zo is 
the unit vector along the z-axis, and I n ,q,m is the ampli- 
tude of the harmonic current I n ,q,m- We are interested 
mostly in the magnetic response of the structure near 
the SRR resonant frequency, ujq = 1/ \JLC§. Our choice 
is motivated by two reasons: (i) nonlinear effects are es- 
sentially enhanced for the frequencies close to the reso- 
nant frequency, and (ii) the left-handed properties of the 
composite of wires and SRRs are observed for the fre- 
quencies on the right-hand side from the SRR resonance. 
Equation Q gives an explicit relation between the am- 
plitudes of the current in SRR, I n ,q,mi and the voltage 
applied across the SRR slit, U n ^ m : 

Z4, g , m = l n , q , m /iu [Co + AC NL (\U g \ 2 )] . (6) 

Next, we introduce the dimcnsionlcss variables r — ujot, 
CI = (u> - Wo)M), and ^ n , q ,m = 1 n ,q,m/1c, where 
X c = looCoUc is the characteristic nonlinear current, 
U c = E c xd g is the characteristic voltage, = Mm ±/L 
are the coupling coefficients, 7 = R/Lloq is the damping 
coefficient, and S n ,q.m = — LoHoira 2 / caJoLZ c is the nor- 
malized electromotive force. The coupling coefficients 
K|i,_L can be calculated numerically for any SRR ge- 
ometry, and the magnetic field of the loop current is 
well known (see, e.g., Ref. However, here we use 

an approximate expression for the magnetic field of the 
loop current |Toj |. which yields k± ss ^(a/d) 3 = k and 
K|| « 2k. In the dimcnsionlcss units, the equations for 
the slowly varying amplitude of the current in SRR (and, 
respectively, the magnetic momenta) can be written as 



i d " - {2n -1-/ + a|*„, 9 , m | 2 } *„, 9 , m - ] 

— 2k {^ n ,q,m + l + ^n,q.m-l — 2ty n , q , m } — K {^' rl +l j 

,q-\-l,m + Vn,q-\,m ~ ^n,q,m} ,(7) 

where t = cuot. 

In addition, Eq. may include long-range effects due 
to the interaction with next-neighboring and other SRRs; 
such effects will produce additional terms in the right- 
hand side of Eq. 10 in the form of higher-order discrete 
derivatives, and they can be shown to remain small not 
changing the major conclusions of our analysis. 

Equation Q has a clear physical meaning. The cou- 
pling coefficient, k, determines a shift of the oscilla- 
tor eigenfrequency due to an effective coupling between 
SRRs, and the nonlinear term produces an eigenfre- 
quency shift due to the nonlinear self-action effect. The 
effect of the higher-order derivatives will be an additional 
shift of the resonant frequency. The right-hand side of 
Eq. J7J represents the second-order difference operators, 
which determine the character of the wave propagation 
and diffraction. The opposite signs between the respec- 
tive discrete difference operators make the derived equa- 
tions fundamentally different from those describing the 
dynamics of two-dimensional discrete systems, where the 
difference operators have the same sign. 
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FIG. 2: (color online) Two-dimensional anisotropic discrete 
diffraction of linear magnctoinductive waves in a lattice of 
SRRs. Shown is the SRR magnetization on the plane (n, m) 
(a) at the upper cutoff frequency when !72 = 3k, and (b) at 
the resonant frequency when Q — 0. 



0.2 




FIG. 3: (Color online) Dependence of the absolute value of 
homogeneous magnetization |$o| of the metamaterial on the 
absolute value of the external field |Eo| for two cases: (a) 
n = 0.2, a = -1, (b) fi = 0.2, a = +1. 



Linear magnetoinductive waves. First, we describe lin- 
ear magnetoinductive waves |]| in lossless metamaterials, 
i.e., we assume d/dr — 0, a = 0, 7 = 0, and X„ jq . m = 0, 
and substitute ^ n ,q,m = F(Cl, k) exp(— ink x d — iqk y d — 
imk z d), where k x ^ y>z are the Cartesian components of 
the wavevector k. The corresponding dispersion relation 
takes the form 



Q = 4k sin 2 



k z d 
2 



- 2k 



sin 



h x d 



and its real solutions describe the linear waves with the 
frequencies |f2| < 4k. For the dimensional frequencies, 
this relatively narrow region is centered at the SRR eigen- 
frequency. We note that the upper cutoff of the frequency 
band of linear waves corresponds to the edge of the Bril- 
louin zone at k x = k y = and k z = n/d, and in this limit 
the linear waves are longitudinal (magnetic plasmons). 
The lower cutoff of the frequency band corresponds to 
the point k z = 0, k x = k y = n/d, and the correspond- 
ing magnetoinductive linear waves are transverse. At the 
resonance (f2 = 0) and in the vicinity of the origin of the 
Brillouin zone, the isofrequency surface forms a resonant 
cone defined as 2(k z d) 2 = (k x d) 2 + (k y d) 2 , and the phase 
velocity of the linear magnetization waves becomes or- 
thogonal to their group velocity. 

To study the generation and propagation of the mag- 
nctoinductive waves in the SRR lattices, we perform nu- 
merical simulations from Eq. J7J). For simplicity, we con- 
sider a two-dimensional geometry, when all fields are ho- 
mogeneous along the y-axis. We simulate a lattice of 
51 x 51 array of SRRs exciting only one SRR in the 
middle of the structure (at n = m = 26), taking the 
damping coefficient 7 = 0.01 and different excitation fre- 
quencies, O = and Q = 3k. To suppress the wave 
reflection from the boundaries of the numerical domain, 
we increase smoothly the damping coefficient in the three 
rows of resonators near the boundaries. At the upper cut- 
off frequency, Q = 3k, the longitudinal wave is shown in 
Fig-EI a )- hi the case of the resonance, SI = 0, we observe 
a resonant cone as shown in Fig. |2fb). 

Nonlinear magnetization. To study nonlinear effects, 
we consider a metamaterial illuminated by a homoge- 
neous electromagnetic wave. The external electromotive 



force is the same in all SRRs: E„ i9)ln = Sq- Steady homo- 
geneous magnetization of the SRR lattice (^ n ,g,m = ^0) 
can be obtained from Eq. (JJJ and is described by the non- 
linear dispersion, which relates magnetization amplitude 
with the external magnetic field: 



{( 



2{l + a\y \ 



\ 2 1 2 
) +7 



}l*o| 



ISol 



(8) 



Two examples of the dependence of vf^ on are shown 
in Fig. |3 When Sla < 0, this dependence becomes multi- 
valued [see Fig. |3ja)] and, as is shown below, the middle 
branch corresponds to unstable solutions. Figure |3|b) 
shows the dependence of \t?o on £0, when Sla > 0. 

Such a bistable dependence of the metamaterial pa- 
rameters may allow switching between the left-handed 
transparent and right-handed opaque states by applying 
a varying external magnetic field and may support 
the existence of dynamically induced transparent regions 
inside an opaque material slab and the formation of spa- 
tiotemporal electromagnetic solitons. 

Modulational instability. Next, we study linear stabil- 
ity of the homogeneous nonlinear states JSJ) with respect 
to small perturbations of the form, 5^„.q.m = (u+v)e xt + 
(u* — v*)e x t , where u = A\ exp(ik x nd + ik y qd + ik z md) 
and v — A2 exp(— ik x nd — ik y qd — ik z md). Substitut- 



ing * r 



*o + &^n, q ,m into Eqs. Q), we obtain the 



following expression for the instability growth rate A, 



A 



7 + J -3|* n | 4 - San\^ \ 2 - 40 2 



(9) 



where £1 = O - 2K[2sin 2 (fc z d/2) - sin 2 (k x d/2) - 
sin 2 (k y d/2)] can take any value in the range [Q — 4k, fl + 
4k] . A simple analysis shows that the real part of A can 
become positive and, therefore, the modulational insta- 
bility can occur when ail < 0. In this case, the bound- 
aries of the modulational instability region are defined 
from the relations. 



|* | 2 = (40/3) ± ^ (4fi 2 /9) - ( 7 2 /3). (10) 

Both the instability and bistability regions can be pre- 
sented on the parameter plane (^0, as shown in Fig.|I| 
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FIG. 4: (Color online) Instability region (bounded by a solid 
curve) on the parameter plane (*l/o,fi)- Dashed: the region 
of the decreasing branch of the homogeneous magnetization 
shown in Fig. 1^1 a). Shaded: instability domain for the magne- 
tization with non-bistable behavior for a£l > [see Fig.|5f b)]. 



ST in such a way that the lower state excites a homoge- 
neous magnetization of the metamaterial corresponding 
to the lower branch of Fig. [21 while the higher state ex- 
cites the upper branch. After a steady state is reached, 
we switch off the inhomogeneous part (put ST = 0), and 
observe the propagation of a nonlinear switching wave 
(kink) with the profile shown in Fig. [SJ Choosing the 
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for the parameters a = — 1 and k = 0.0025. The case 
a = +1 can be recovered from the same results by chang- 
ing f2 — > — O. The curve that depicts the instability re- 
gion consists of three parts. The upper part represents 
the boundary (|10|) for Cl = SI + 4k, which corresponds 
to the instability with respect to the excitation of the 
transverse waves (at k z = 0, k x = k y = n/d). The lower 
curve corresponds to the excitation of the longitudinal 
waves (at k z = it/d, k x = k y = 0). Along the horizon- 
tal line the instability occurs for the waves propagating 
at some angle to the main axis of the lattice. Thus, the 
instability covers completely the region corresponding to 
the middle branch of the bistable dependence shown in 
Fig-EI a )- Some parts of the upper and lower branches of 
that dependence are also unstable. The instability also 
occurs for aQ, > 0, as shown in Fig. 0] (shaded region). 

Existence of unstable regions can be crucial for engi- 
neering tunable mctamatcrials with nonlinear properties. 
The magnetization waves propagating inside the meta- 
material represent alternating positively and negatively 
magnetized SRRs. As a result, the average magnetiza- 
tion of the metamaterial will be dramatically reduced, 
possibly suppressing the regions of the negative effective 
permeability. 

Domain walls. Finally, we study numerically nonlin- 
ear switching waves (or domain walls) which can prop- 
agate in the nonlinear metamaterial. We solve Eq. Q) 
numerically for a two-dimensional SRR lattice with the 
size NxN, when N = 100, for defocusing nonlinearity 
(a = —I) and = 0.2, when no linear waves exist. First, 
we create an inhomogeneous magnetization by a step-like 
external field T n = T + ST x x(n — N/2), where x = — 1 
for n < N/2, and x = +lj otherwise. We chose Eq and 



FIG. 5: (color online) Example of a switching wave of mag- 
netization (kink) in a nonlinear metamaterial. Shown is the 
magnetization = 'I'n.o.o- Arrows show schematically the 
magnetization of the corresponding parts of the metamaterial, 
i.e., the magnetization is positive for arrows directed upwards, 
and negative, otherwise. 



parameters Tq and ST, we may control the speed and di- 
rection of the kink propagation. Two out-of-phase kinks 
may create a magnetization domain which differs from 
the rest of the material. Depending on the external field, 
such a domain can collapse, expand, or preserve its shape. 
The possibility to control creation and dynamics of such 
domains seems to be promising for the design of the struc- 
tures with controllable periodic magnetization, photonic 
crystals, which parameters can be made tunable. 

In conclusion, we have analyzed, for the first time 
to our knowledge, nonlinear magnctoinductive waves in 
composite left-handed mctamatcrials. We have derived 
the effective discrete model that describes the propaga- 
tion of magnetoinductive waves in a lattice of split-ring 
resonators, and studied both linear and nonlinear waves, 
demonstrating they can change dramatically the aver- 
age magnetization of left-handed materials. We have re- 
vealed that the bistable magnetic response may lead to 
the propagation of domain walls in nonlinear mctamatc- 
rials. 
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